Ultra-fast propagation of Schrodinger waves in absorbing media 
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We identify the characteristic times of the evolution of a quantum wave generated by a point source 
with a sharp onset in an absorbing medium. The "traversal" or "Biittiker-Landauer" time (which 
grows linearly with the distance to the source) for the Hermitian, non-absorbing case is substituted 
by three different characteristic quantities. One of them describes the arrival of a maximum of 
the density calculated with respect to position, but the maximum with respect to time for a given 
position becomes independent of the distance to the source and is given by the particle's "survival 
time" in the medium. This later effect, unlike the Hartman effect, occurs for injection frequencies 
under or above the cut-off, and for arbitrarily large distances. A possible physical realization is 
proposed by illuminating a two-level atom with a detuned laser. 
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As it was well understood by Brillouin and Sommer- 
feld long ago, certain wave features may travel at ve- 
locities exceeding c in systems described by relativistic 
equations, without violating Einstein's causality 0. In 
non-relativistic equations one may similarly find "ultra- 
fast" phenomena which are subject to non-relativistic 
causality conditions In particular, the Hartman effect 
has been studied thoroughly theoretically and exper- 
imentally both for non-relativistic and relativistic equa- 
tions H S H IS II H H ES EH : the peak of the transmit- 
ted packet in a collision of a particle with a square barrier 
emerges from the barrier edge at a time independent of 
the barrier width d when the incident (average) energy is 
below the barrier energy, and provided that the barrier is 
neither too thin nor too opaque. This later large-d lim- 
itation occurs because, for very opaque conditions, the 
wave packet components above the barrier dominate and 
lead to a linear-with-c? behaviour of the time of the peak 
0, For the barrier lengths where the Hartman ef- 
fect holds one could in principle transmit a peak in a 
waveguide in such a way that it arrives simultaneously 
to receivers located at different distances. We describe 
in this letter an even stronger effect for the Schrodinger 
equation: if the particles are sent by a point source with 
a sharp onset, or by a sudden opening of a shutter, to a 
region with an effective absorbing potential, the tempo- 
ral maximum arrives at equal times at arbitrarily large 
positions, i.e., without the limitation of the Hartman ef- 
fect to moderately opaque conditions. Also at variance 
with the Hartman effect, this ultra-fast propagation of 
a wave feature due to absorption also holds for injection 
energies above the potential threshold. 

Consider, for x > 0, the following Schrodinger equa- 
tion with "source" boundary conditions (all quantities 
are dimensionless unless stated otherwise), 



.dtp 



d 2 ip 
dx 2 



+ (1 - iVM 



at 

i/>(x = Q,t) =e- jMot ®(t) 
ip(x,t)=0 x>0,t<0. 



(1) 



where ujq is the injection frequency. This is a non- 
Hermitian generalization of the source with a sharp onset 
studied by several authors before |H 111 111 111 111 Il7| . 
The imaginary potential —iV\ models the passage from 
the incident channel to other channels which are not rep- 
resented explicitly (a physical example is provided be- 
low). The dispersion relation corresponding to Eq. Q 
is 



k = y/u- l+iVx, Im(fc) > 0, 
and the solution to Eq. (JTJ is given by 

Akx—iut 



ip(x,t) = 
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duj — -, 

^ LJ LJq + 7,0 



or, in the fe-complex plane, as 



where 
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fc(wo) = v^o- 1 + iVi, Im(fco) > 0, (6) 



and the countour Tk goes from — oo to oo passing above 
all singularities. Introducing a new integration variable 



1 + i 
V2 



(7) 



Eq. Q takes the form 

i e ~iVt+ix 2 /it 

tp(x,t) = - 
where 
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FIG. 1: \ip\ versus x for three different values of t: t — 250 
(solid line), t = 300 (dashed line) and t — 500 (long dashed 
line). We have chosen uio = 0.99 and Vi = 0.001. 



FIG. 2: \ip\ versus time for three different values of V\\ Vi = 
(solid line), Vi = 0.001 (dashed line) and Vi = 0.01 (long 
dashed line). The position has been fixed to x = 100 and 
lo = 0.9. 



and the contour r„ goes from — oo to oo passing above 
the two simplepoles. Using the integral definition of the 
w— functions [lfll^j Eq. © is finally given by 



2 /4t 



[w(-uq) + w(-u' )] . 



(10) 



To interpret this exact solution in simple terms it is 
useful to find approximations. In particular, the contour 
may be deformed along the steepest descent path from 
the saddle at x/2t in the fc-plane (the origin in the u- 
plane). The path crosses the pole at kg at a time 



a;/2[Re(fc )+Im(fco)]- 



(11) 



Eq. JSJ can be written as the sum of contributions from 
the saddle point and the pole, 



ip p (x,t) 
ijj s (x,t) 



4> P (x,t) +ip a (x,t) 
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(l+i)k (r 2 +t 2 ) 



where 
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(12) 
(13) 

(14) 
(15) 

(16) 



Alternatively, Eq. (|12fl may be obtained from the asymp- 
totic expansion of the w's for large values of the mod- 
ulus of their arguments \u \ and \u' \. They become 



large for large |fcrj|, x or t, and also for very small t. 
Both have minima at t = |r|, with minimum values 
{x[\k \ ± Re(fc )]} 1/2 - (A scale for the validity of the 
saddle plus pole approximation for all t is thus x > 



[|fco| — Re(ko))" 1 ). For injection frequencies below the 
cut-off the saddle contribution dominates up to exponen- 
tially large times as in the Hermitian case [13, so that 
r c is not of much significance whereas, above threshold, 
t c is a good scale for the arrival of the main signal. The 
discussion hereafter refers to the case luq < 1 (injection 
below cut-off) unless stated otherwise. 

In Figure n we can see the formation of the forerunner 
in a sequence of three snapshots of the density at three 
different instants and for a fixed value of V±. (For in- 
creasing values of Vj the spatial peak appears at larger 
x and it also takes a longer time to be formed.) From 
d\ip(t, x)\ 2 /dx = we may obtain the position x(t) of 
the "spatial" maximum at time t. This also defines (by 
inverting x(t)) a function ts(x), namely, the time when 
this spatial maximum arrives at x. One finds from Eq. 
(|14fl that the spatial maximum in the large- a; region is 
given by 



ts 



(17) 



a role played by the real r in the absence of absorption 

E3. 

A single real quantity t(V\ = 0) in the Hermitian case 
has been substituted, for a non-zero Vi , by three different 
quantities: the time for pole-cutting t c , a complex r, 
and its modulus |r|, all of which tend to the "Biittiker- 
Landauer" traversal time r(V"i = 0) without absorption 

Similarly, we may fix x, calculate d\ip(x, t)\ 2 /dt = 0, 
and solve for t to obtain a "temporal" maximum, tt(x). 

In Figure |21 we have plotted the density, \ip\ 2 , versus 
time for a large x, so that, in the scale used, the exact 
solution and the saddle aproximation ip s are undistin- 
guishable. One clear effect of the complex potential is 
the decrease of the amplitude; also, the peak arrives ear- 
lier when V\ increases. At variance with the Hermitian 
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case, the time of arrival of the maximum is not propor- 
tional to t (jt = t/V3 when V\ — 0). The equation 
d\ij) s \ 2 /dt = cannot be solved analytically for V\ ^ 
in a generic case, so there is no explicit formula for vp. 
Nevertheless, if x 3> 2 1 / 2 |fc |/Vi, one finds 



tt 



1 



(18) 



i.e., the temporal maximum coincides with the mean sur- 
vival time of a particle immersed in the absorbing poten- 
tial, and it is independent of x and loq, which is the most 
important result of this work. A time- frequency analysis 
|l6l |21| shows that the maximum is not tunnelling but 
it is dominated by frequencies above the cut-off, in par- 
ticular by the frequency 1 + (xVi) 2 corresponding to the 
"classical" velocity required to arrive at x at time tt{x). 
In fact, this effect is also present for injection frequencies 
above the cut-off. It has, nevertheless, no purely classical 
explanation, in the sense that any classical ensemble of 
particles injected at a constant rate into the absorbing 
medium with an arbitrary momentum distribution from 
t = on, would lead, for fixed x, to a monotonous in- 
crease of the density up to the asymptotic, stationary 
value. 

The important difference between temporal and spatial 
maxima may be understood from the scaling satisfied by 
the saddle term, 



e -2tvi(»?-i) 
\^ s (f]x lV t)\ 2 = |Va(M)| S 



(19) 



Since the exponential depends on t, but not on x, the 
spatial maximum travels at constant velocity whereas, 
the temporal maximum does not, except for V\ = 0. 

Some recent works have examined the time of arrival of 
the maximum at a given position for the source problem 
without absorption also for small values of x, where the 
pole-saddle approximation is not valid |17ll22j . For small 
x, tt{x) presents in that case a basin with a minimum. In 
Fig. we have plotted this quantity versus x for differ- 
ents values of the absorbing potential V\ . Note that the 
basin dissappears by increasing the absorption. The most 
prominent feature of the curves though is their constant 
value 1/2 Vi for large enough x, instead of the linear de- 
pendence found without absorption. Unlike the Hartman 
effect, the arrival of the temporal peak in an absorbing 
medium stays constant for arbitrarily large x. The effect 
is also present for other boundary conditions, in particu- 
lar for the "Moshinsky shutter" boundary condition p3| 
corresponding to the initial, truncated-plane-wave state 



1 



tp(x,t = 0) = -— e lkox Q(-x) 



2tt 



and the potential 



V(x) = (l-iVx)Q(x). 



(20) 



(21) 



Using the same techniques applied in [24(, we may cal- 
culate ip(x,t). The corresponding density is shown in 
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FIG. 3: (a) Time of arrival of the temporal peak versus x for 
Vi = (solid line), V\ = 0.001 (dashed line), Vi = 0.01 (long 
dashed line) and V\ =0.1 (dotted-dashed line). The injection 
frequency is luq = 0.81. (b) Time of arrival of the peak versus 
x for Vi = (solid line), Vi = 0.01 (dashed line), Vi =0.1 
(long dashed line) and Vi = 0.5 (dots), and ko = 0.9. The 
circle, rhomb, and triangle mark the value of 2 1 / 2 |feo|/Vi for 
Vi = 0.1,0.01, and 0.001, respectively. 



Fig. |2t>, where the same characteristic time 1/2 Vi which 
describes the arrival of the peak at large x is found. 

A close physical realization of the Moshinki shutter 
problem with absorption may be based on a fluoresecence 
experiment where an atom is first prepared according to 
a truncated plane wave in an internal state |1), and then 
let evolve after a sudden shutter opening at x = t = 0. 
The "external" region, x > 0, is illuminated with a per- 
pendicular laser. Let us assume a A— configuration for 
three relevant atomic levels such that the laser couples 
levels |1) and |2) whereas |2) decays irreversibly by spon- 
taneous photon emission to a ground state |0) . According 
to the quantum jump technique |25j the amplitudes for 
levels |1) and |2) obey an effective Schrodinger equation 
with Hamiltonian (all quantities are now dimensional) 



H = p 2 /2m 




o n@(x) 



(22) 



where f2 is the Rabi frequency, A the detuning between 
the laser frequency and the transition frequency wi 2 , 7 is 
the inverse life time of |2), and p the momentum operator 
for the initial atomic direction X. For large detuning, 
|2A + i7| >> ft, a further reduction is possible to an even 
simpler effective Schrodinger equation for the amplitude 
of level 1 with potential |2!J, yCJ, |31| 



V{X) = 



hAn 2 o(x) .h 1 n 2 Q{x)/2 

4A 2 + 7 2 



(23) 



4A 2 + 7 2 ' 

which can be transformed by appropriate scaling to the 
form of Eq. (|21fl . We have solved numerically the full 



4 



two-channel Schrodinger equation with the Hamiltonian 
of Eq. and the initial state of Eq. (|2o"|) for the 

1 1)— component. We have found tt(x) = l/2Vi for ar- 
bitrarily large distances to the shutter, in full agreement 
with the solution of the one-channel equation with ab- 
sorbing potential, Eq. QJ. 
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